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SECTION  1 
INTRODUCTION 


Over  the  past  several  years  the  interactions  between  spacecraft  and  the  vacuum, 
electromagnetic  and  particle  radiation  environment  in  space  have  been  of  considerable 
interest  to  the  scientific  and  engineering  community.  Electrostatic  charging  of  high 
altitude  satellites  by  the  effects  of  the  ambient  plasma  and  solar  radiation  is  of  serious 
importance,  since  it  can  affect  the  operation  of  several  on-board  systems  in  scientific, 
military,  and  commercial  satellites. 

Some  laboratory  facilities  are  presently  available  (see  Refs.  1-3,  for  example)  and 
(4) 

others  are  projected  which  can  provide  a  degree  of  simulation  of  the  radiation  and 
plasma  environment  of  space.  A  simulation  facility  is  necessarily  limited  by  engineering 
and  economic  considerations  to  providing  only  a  partial  re-creation  of  the  actual 
environment  encountered  by  a  space  vehicle.  For  example,  only  a  few  monoenergetic 
beams,  incident  on  a  target  from  a  limited  number  of  directions,  can  be  envisioned  to 
simulate  geomagnetic  substorm  plasmas  which  generally  have  distributed  energy  spectra 
and  are  fairly  isotropic.  The  question  as  to  what  constitutes  a  "good"  simulation  of  a 
given  plasma  environment  has  often  been  posed,  but  there  are  presently  no  consistent 
answers. 

The  purpose  of  the  present  investigation  is  to  analyze  specific  techniques  for 
choosing  the  energies  and  current  densities  of  monoenergetic  electron  and  ion  beams  to 
simulate  given  space  plasmas  with  distributed  particle  energy  spectra.  The  plasmas  are 
characterized  by  various  averages  of  their  velocity  distribution  functions,  and  the 
energies  and  particle  densities  of  the  monoenergetic  beams  are  computed  so  that  there  is 
c  good  match  between  the  characteristics  of  the  distributed  spectra  and  those  of  the 
beams.  We  hope  that  this  analysis  will  help  provide  a  concrete  link  between  the 
measurements  of  the  space  environment  and  the  engineering  specifications  for  charged 
particle  beams  in  laboratory  simulation  facilities. 

This  interim  report  outlines  the  approach  which  we  have  followed  to  define  beam 
energies  and  densities  to  simulate  distributed  plasma  spectra.  It  does  not  consider  the 
equally  important  problems  of  the  effects  of  discrete  angular  directions  of  the  charged 
particle  beams  or  the  interactions  of  the  incident  beams  with  the  target.  We  have, 
however,  made  considerable  progress  in  both  these  areas,  which  will  be  reported  in  other 
publications  and  in  the  final  report  of  the  program. 
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2.1  BACKGROUND 

A  plasma  can  be  characterized  by  various  averages  of  the  velocity  distributions  of 
its  constituent  particles.  In  general,  the  "velocity  moments"  of  a  given  distributon 
function,  f(v),  are  defined  by 

oo 

J'  f(v)  v^dv  ^ 

k  =  0,1,2, - . 

o 

where  the  4  jrv  dv  term  represents  an  infinitesimal  element  in  (isotropic)  velocity  space. 

The  velocity  moments,  M^,  can  be  related  to  physical  averages  for  several  values 
of  k.  For  example,  MQ,  M^,  M2,  and  M3  are  related,  respectively,  to  the  average  number 
density  <N>,  particle  flux,  <NF>,  pressure,  <P>,  and  energy  flux,  <EF>,  of  the  given 


jj 
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particle  type  in  the  plasma. 


The  average  speed,  <v>,  in  Equation  (3)  is  defined  by 


<v> 


(8) 


The  expressions  on  the  right-hand  side  of  Equations  (2) -(7)  are  given  for  the  case  of  a 
Maxwellian  velocity  distribution, 


f(v) 


(9) 


where  n,  m,  and  T  are  respectively  the  number  density,  mass,  and  temperature  of  the 
particles  and  k  is  Boltzmann's  contant. 

Average  and  RMS  "Temperatures11 

A  useful  method  for  characterizing  a  non-Maxwellian  plasma  is  to  define  effective 

(5) 

temperatures  which  are  related  to  ratios  of  the  velocity  moments.  The  average  and 
RMS  temperatures  are  given  by 


T  1  -  m  M2 

AV  k  <N>  3k  MQ 


(10) 


=  i  <EF>  =  m_  “3 
k  2<NF>  4  k  Mj 


(11) 


The  two  temperatures  are  equal  when  the  velocity  distribution  is  Maxwellian. 

2.2  MONOENERGETIC  BEAMS  TO  MATCH  VELOCITY  MOMENTS 

A  technique  to  simulate  a  plasma  with  a  distributed  velocity  distribution  is  to 
choose  the  velocities  and  particle  densities  of  mononergetic  beams  so  that  their  velocity 
moments  match  those  of  the  plasma.  Under  these  conditions,  the  average  parameters  of 
the  beams,  such  as  number  density,  pressure,  or  energy  flux,  are  equal  to  those  of  the 
plasma  component  under  simulation. 


In  general,  a  single  beam  can  match  two  moments  of  the  distributed  spectrum,  so 
that  two  beams  can  match  four  moments,  three  beams,  six  moments,  etc.  As  discussed 
in  Section  2.4,  it  is  also  possible  to  overspecify  the  problem  and  use  more  than  the 
minimum  number  of  beams  to  match  a  given  number  of  velocity  moments. 

2.2.1  Single  Beam  Energy 

A  monoenergetie  beam  can  match  two  moments  according  to  the  simultaneous 
equations, 


n,  v, -1  =  M 
b  b  3 


k  _  u 

nbVb  -  “k 


(12) 


where  nb  and  v^  are  the  density  and  velocity  of  the  beam  particles. 

For  example,  when  the  zeroth  (number  density)  and  second  (pressure)  moments  are 
chosen, 


(13) 


or,  in  terms,  of  beam  energy,  E^, 


If  the  first  (number  flux)  and  third  (energy  flux)  moments  are  used, 


(14) 


(15) 


-4- 


or 


E,  =  2  k  T 


RMS 


(16) 


2.2.2  Two  Beam  Energies 

The  densities  and  velocities  of  two  monoenergetic  beams  can  be  found  to  match 
the  zeroth  through  third  velocity  moments  of  the  distributed  spectrum  by  solving  four 
simultaneous  equations: 


nl  +  n2  =  n 


nlVl  +  n2V2  =  n  <V> 


"lVl2  +  n2V22  =  I  n  TAV 


nlVl3  +  n2V23  =  i  n  <V>  TRMS 


(17) 


where  n^  r\y  Vj,  and  Vg  are  the  densities  and  velocities  of  the  beams,  and  the  velocity 
moments  of  the  distributed  spectrum  have  been  replaced  by  relations  (8),  (10),  and  (11). 
Boltzmann's  constant,  k,  has  been  taken  to  be  unity. 

It  is  shown  in  Appendix  A  that  the  velocities  and  densities  of  the  monoenergetic 
beams  can  be  found  analytically. 
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(19) 


For  a  Maxwellian  plasma,  where  T ^  =  TRMg  =  T,  Eq.  (18)  simplifies  somewhat, 


The  beam  densities  and  energies  are  then  found  to  be 

n  =  0.382  n 

(21) 

=  0.618  n 


c1  =  3.007  T 
E2  =  0.568  T 


(22) 


2.2.3  Three  Beam  Energies 

Six  moments  of  the  distributed  spectrum  can  be  used  to  compute  the  densities  and 
velocities  of  three  monoenergetic  beams.  No  analytical  solutions  have  been  found  for 
this  case,  but  iterative  techniques  can  be  used  to  find  solutions  of  the  set  of  six 
simultaneous,  nonlinear  equations. 

As  discussed  in  Appendix  B,  the  beam  velocities  and  densities  can  be  found  in 
terms  of  the  average  speed  and  density  of  the  plasma  particles.  For  the  case  of  a 
Maxwellian  plasma  with  temperature,  T,  the  beam  densities  and  energies  are 


n  =  [0.087,  0.588,  0.325]  n 

E:  2  3  =  [6.931,  1.657,  0.303]  T 


(23) 
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Different  values  will  be  found  for  other  types  of  velocity  distribution  functions, 
but  the  method  used  to  compute  the  Maxwellian  results  is  general  for  all  realistic 
spectral  shapes. 

2.3  TWO-MAXWELLIAN  PLASMAS 

Garrett  showed  that  a  two-Maxwellian  fit  is  often  a  good  representation  of  plasma 

(61 

distribution  functions  measured  during  geomagnetic  substorms.'  7  The  density  and 
temperature  of  each  Maxwellian  component  can  be  found  from  four  velocity  moments  of 
the  measured  spectrum.  It  is  possible,  in  principle,  to  find  three-Maxwellian  fits  which 
match  six  moments,  although  the  effects  of  errors  in  measurement  of  the  plasma 
spectrum  become  increasingly  exaggerated  when  computing  the  high-order  moments.  It 
should  also  be  possible  to  find  multiple-Maxwellian  least-square  fits  directly  from  the 
measured  distribution  functions  without  computing  the  velocity  moments  of  the  data. 

2.3.1  Single  Beam  Energy 

A  two-Maxwellian  distribution  has  average  and  RMS  temperatures  given  by 
_  nlTl  +  n2T2 

fAV  nx  +  n2  (24) 


3/2  3/2 

W  =  1/2  "  1/2 

nlTl  +  n2  z 


(25) 


where  n^,  n2,  Tj,  and  T2  are  the  respective  densities  and  temperatures  of  the  two 
components  of  the  spectrum. 

A  single  monoenergetic  beam  can  match  two  velocity  moments  of  the  distributed 
spectrum  if  its  density  and  energy  are  chosen  according  to  Eqs.  (13) -(16)  above.  For 
example,  if  the  beam  density  is  equal  to  the  total  plasma  density,  nj  +  n2,  and  its  energy 
is  3/2  TAy,  then  the  zeroth  and  second  velocity  moments  of  the  two-Maxwellian  plasma 
and  the  monoenergetic  beam  are  equal. 


-7- 


2.3.2  Two  Beam  Energies 


Two  methods  exist  for  matching  the  velocity  moments  of  a  two-Maxwellian 
distribution  by  two  monoenergetic  beams.  First,  the  energy  and  density  of  each  beam 
can  be  chosen  individually  to  match  two  moments  of  each  of  the  Maxwellian  components 
of  the  spectrum.  In  this  case,  Eqs.  (13) -(16)  would  be  employed  along  with  the  densities 
and  temperatures  of  the  two-Maxwellian  fit. 

The  second  approach  is  to  use  the  average  and  RMS  temperatures  of  the 
two-Maxwellian  fit,  Eqs.  (24)  and  (25)  and  to  calculate  the  beam  velocities  and  densities 
from  Eqs.  (18)  and  (19).  In  both  cases,  as  many  as  four  moments  of  the  two-Maxwellian 
distribution  function  can  be  matched  by  two  monoenergetic  beams.  In  practical 
situations,  physical  considerations  would  be  required  to  make  a  choice  between  the  two 
methods  of  matching  velocity  moments. 

2.3.3  Three  or  More  Beam  Energies 

The  moments  of  a  two-Maxwellian  distribution  function  can  be  matched  in  several 
different  combinations  with  multiple  monoenergetic  beams.  As  in  the  two-beam  case, 
each  Maxwellian  component  of  the  plasma  can  have  one  or  more  beams  assigned  to  it 
which  individually  match  velocity  moments.  For  six-moment  matching,  three  beam 
energies  and  densities  could  be  selected  using  Eq.  (23)  for  each  component,  and  a  total  of 
six  beam  energies  would  be  required  to  simulate  the  two-Maxwellian  plasma.  As 
mentioned  in  Section  2.2.3,  the  computed  values  of  the  zeroth  through  fifth  moment  of 
the  full  spectrum  can  also  be  used  directly  to  find  three  beam  energies  and  densities 
through  the  iterative  minimization  procedure  described  in  Appendix  B. 


2.4  ARBITRARILY  ASSIGNED  BEAM  ENERGIES 


The  velocity  moments  of  a  measured  distribution  function  can  also  be  matched  by 
monoenergetic  beams  whose  velocities  are  chosen  arbitrarily.  As  an  example,  four 
monoenergetic  beams  can  match  four  velocity  moments: 


nl  +  n2  +  n3  +  n4  =  M0 


v.n,  4-  v„n„  +  v_n_  +  v.n,  =  M, 
11  22  33  44  1 


Vl2nl  +  V22n2  +  V32n3  +  V42n4  "  M2 


(26) 


3  3  3  3 

V1  nl  +  v2  n2  +  v3  n3  +  v4  n4  =  M3 


When  the  four  beam  velocities  are  fixed,  then  it  is  only  a  matter  of  solving  a  set 
of  linear  simultaneous  equations  for  the  beam  densities,  n1  through  n4«  It  should  be 
pointed  out  that  not  all  combinations  of  beam  velocity  may  be  chosen  for  Vj  through  v^, 
because  negative,  and  therefore  unphysical,  solutions  for  the  beam  densities  can  be 
obtained  in  some  cases. 

Table  1  gives  the  densities  calculated  for  three,  four,  and  five  beams  as  a  function 
of  preassigned  beam  energies.  The  beam  energies  and  densities  are  normalized  to  the 
temperature  and  density  of  a  Maxwellian  distribution,  and  the  velocity  moments  used  for 
the  calculations  are  given  by  the  right-hand  side  of  Eqs.  (2) -(6).  The  first  three-beam 
solution  in  Table  1  is  a  check  of  the  six-moment  solution,  Eq.  (23),  found  by  the  iterative 
procedure  discussed  in  Section  2.2.3. 
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TABLE  1.  PARTICLE  DENSITIES  FOR  PREASSIGNED 
BEAM  ENERGIES 


SECTION  3 

CONCLUSIONS  AND  RECOMMENDATIONS 


It  has  been  shown  in  Section  2  that  the  low  order  velocity  moments  of  plasmas 
with  arbitrary  distributed  energy  spectra  can  be  matched  by  combinations  of 
monoenergetic  beams.  The  beams  simulate,  in  a  mathematically  precise  manner,  certain 
properties  of  an  undisturbed  space  plasma.  However,  interactions  between  the  plasma 
and  a  complicated  object  such  as  a  satellite  were  not  considered  in  this  analysis. 
Spacecraft  charging,  on  the  other  hand,  is  a  phenomenon  which  depends  not  only  on  the 
properties  of  the  undisturbed  plasma,  but  also  on  the  detailed  interactions  of  the  plasma 
particles  with  the  surface  and  the  magnitude  and  shape  of  the  equipotentials  around  the 
spacecraft. 

(7) 

Charging  models  of  various  degrees  of  complexity  have  been  developed  in  an 
attempt  to  account  for  the  observed  interactions  between  the  environment  of  space  and 
instrumented  spacecraft.  An  important  task  of  the  present  program  is  to  use  such  a 
model  to  study  the  consequences  of  substituting  monoenergetic  beams  for  plasmas  with 
distributed  spectra.  We  believe  that  comparisons  of  this  kind  will  provide  insight  into  the 
selection  of  a  few  combinations  of  beam  parameters  which  will  provide  good  simulation 
of  spacecraft  charging  phenomena  among  the  many,  mathematically  equivalent, 
techniques  for  matching  averaged  properties  of  plasmas  and  monoenergetic  beams. 
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APPENDIX  A 

TWO-BEAM  SOLUTION  TO  MATCH  FOUR  VELOCITY  MOMENTS 
(6) 

Garrett  '  has  shown  for  a  similar  problem  that  the  set  of  nonlinear  simultaneous 
equations,  Eq.  (17),  can  be  reduced  to  a  single  quadratic  equation  of  the  form, 

AX2  +  BX  +  C  =  0  (A-l) 


where  A,  B,  and  C  are  functions  of  the  right-hand  side  of  Eqs.  (17).  For  the  problem  of 
calculating  the  velocities  and  densities  of  monoenergetic  beams,  we  have  set 


(A-2) 


The  beam  velocities  are  found  to  be 


<v> 


1,2 


6  T,„  -2m  <v> 
AV 


2  I4  TRMS  "  3  TAV 


8  trms  (2  trms  "  9  tav  +  2  m  <v>  ) 


(A-3) 


-  27  tav  (‘ - 
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where  Vj  (vj)  corresponds  to  the  +  (-)  sign  of  Eq.  (A-3). 
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The  beam  densities  corresponding  to  the  velocities  are  found  by  substitution, 


(A-4) 
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APPENDIX  B 

THREE-BEAM  SOLUTION  TO  MATCH  SIX  VELOCITY  MOMENTS 


To  find  a  three-beam  solution  to  match  six  velocity  moments  of  a  distributed 
spectrum,  six  nonlinear  simultaneous  equations  must  be  solved.  These  equations  are  of 
the  form 


3 
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k=  1 
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£ 

j=0 


nkvk 


M. 

J 


(B— 1 ) 


velocity  moment  of  the  distributed  spectrum. 

Although  any  physically  reasonable  distribution  function  can  be  utilized,  we  used 
the  six  velocity  moments  of  a  Maxwellian  distribution,  Eqs.  (2) -(7),  to  find  solutions  of 
Eq.  (B-l).  The  beam  densities  and  velocities  were  normalized  by  changing  variables, 


a,b,c 


1.2,3 


n 


x,y,z 


V 1,2,3 
<v> 


so  that  the  moment  equations  for  a  Maxwellian  distribution  become 


(B— 2) 


a  +  b  +  c  =  1 
ax  +  by  +  cz  *  1 

ax ^  +  by^  +  c 2}  “  (B-3) 

ax^  +  by'*  +  cz"*  =  j 

4.4  A  15  2 

ax  +  by  +  cz  =  n 

5  .  .5  .  5  32 

ax  +  by  +  cz  ■  g  it 
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Eqs.  (B-3)  can  be  solved  by  an  iterative  procedure  that  minimizes  the  expression 

represented  by  the  absolute  sum  of  Eqs.  B-3.  Trial  solutions  are  substituted  into  a 

computer  program  which  then  converges  on  the  best  solution  through  an  iterative 

(8) 

process. 

The  only  solutions  of  Eqs.  (B-3)  found  by  the  iterative  routine  were  permutations 
of  the  following: 


(x,y , z)  =  (1.9679,  1.14085,  0.4881) 
(a, b, c)  =  (0.0866,  0.5879,  0.3255) 


(B-4) 


The  energies  of  the  beams  were  found  from 


E 


1 


2 

x 


and  similar  expressions  for  E 2  and  Eg 


(B-5) 
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